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LONG-TERM  GOALS 

My  long-term  goals  are  to  improve  our  understanding  of  the  physical  evolution  of  the  continental  shelf 
and  slope,  and  to  enhance  our  ability  to  extract  stratigraphic  information  from  geophysical  data  of 
these  regions. 

OBJECTIVES 

The  specific  objectives  of  this  project  are  to: 

(1)  Aid  the  development  of  the  integrative  stratigraphic  modeling  system  SEQUENCE  (Steckler,  1999) 
by: 

a)  helping  create  a  component  model  for  the  long-term  evolution  of  shelf  strata,  and 

b)  creating  a  component  model  for  the  long-term  evolution  of  slope  strata. 

(2)  Develop  a  model  for  generating  synthetic  seismograms  of  stratigraphic  simulations  produced  by 
SEQUENCE  and  by  SEDFLUX  (Syvitski  et  ah,  1999),  the  second  integrative  stratigraphic 
modeling  system  being  developed  for  ONR. 

(3)  Continue  work  on  modeling  of  unconfined  submarine  debris  flows  and  validate  this  work  against 
experimental  flows. 

APPROACH 

SEQUENCE  is  being  designed  to  rapidly  simulate  the  first-order  evolution  of  continental  margin 
morphology  over  long  geologic  time  periods.  To  meet  this  requirement,  the  shelf  and  slope  component 
models  for  SEQUENCE  are  1-D,  Eulerian,  finite-difference  models  that  produce  2-D,  dip-line 
simulations  of  shelf  and  slope  evolution.  Both  are  relatively  fast  and  encompass  the  physiographic 
provinces  in  their  entirety,  using  annual  to  multi-year  time  steps,  and  distance  steps  of  tens  of  meters. 
The  shelf  model  simulates  the  formation  of  shelf  strata  under  the  combined  effects  of  sediment  supply, 
waves,  and  a  long-term,  net-offshore  current.  The  slope  model  simulates  episodic  failure  and/or  the 
bypassing  of  sediments  to  the  deep  sea. 
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The  modeling  of  synthetic  seismograms  from  strata  generated  by  SEQUENCE  and  SEDFLUX  requires 
two  separate  models.  The  first  model  generates  the  physical  properties  needed  to  produce  synthetic 
seismograms,  principally  bulk  and  shear  moduli.  The  second  model  uses  these  properties  to  simulate 
the  theoretical  response  of  the  computer-generated  strata  to  seismic  energy,  including  its  attenuation. 

Submarine  debris  flows  and  the  deposits  they  produce  are  being  modeled  in  three  dimensions  using  a 
2-D,  layer-averaged,  numerical  model  of  a  Bingham  fluid.  The  model  is  to  be  validated  against 
unconfined  experimental  debris  flows,  the  heights  and  velocities  of  which  have  been  well  documented. 

The  shelf  and  slope  modeling  is  being  done  in  collaboration  with  Drs.  Pat  Wiberg  (UVA),  John 
Swenson  (UMinn-Duluth),  Chris  Paola  and  Gary  Parker  (UMinn-Minneapolis).  The  seismic  modeling 
is  being  done  with  Dr.  Anastasia  Stroujkova  (Duke/UConn).  And  the  debris  flow  modeling  is  being 
done  with  Drs.  Jim  Buttles  (MIT)  and  Gary  Parker  (UMinn-Minneapolis). 

WORK  COMPLETED 

The  shelf  model  continues  to  be  refined,  but  a  working  version  is  up  and  running.  In  the  model, 
sediments  introduced  at  the  coast  are  kept  in  suspension  by  wave  energy  and  gradually  moved 
basinward  until  water  depths  become  deep  enough  that  the  wave  energy  is  too  low  to  prevent  long¬ 
term  sediment  accumulation. 

The  slope  model  is  also  up  and  running.  In  this  model,  sediment  that  escapes  the  shelf  accumulates  on 
the  slope  until  the  slope  reaches  a  critical  angle.  At  this  point,  the  sediments  either  bypass  the  slope  or 
the  slope  fails  depending  on  the  user's  preference.  The  bypassed  or  remobilized  sediments  are  then 
transferred  to  the  deep  sea. 

In  order  to  test  the  performance  of  the  shelf  and  slope  models,  the  two  have  been  linked  with 
component  models  for  the  coastal  plain  and  the  deep  sea.  The  coastal  plain  model,  which  uses 
diffusion  to  approximate  fluvial  deposition  and  erosion  (Paola  et  al.,  1992),  provides  the  sediment 
input  to  the  shelf  model.  The  deep-sea  model  uses  a  diffusion  equation  derived  from  the  conservation 
equations  for  a  turbidity  current  (Parker  et  al.,  1986)  to  distribute  the  sediment  received  from  the  slope 
model  across  an  adjoining  continental  rise. 

The  physical  property  model  needed  for  generating  synthetic  seismograms  is  complete  (Herrick,  2001; 
Stroujkova  et  al.,  in  prep.).  It  can  be  used  to  simulate  bulk  and  shear  moduli  not  only  for  computer 
simulations  of  strata,  but  also  for  experimental  strata  created  in  the  laboratory.  The  model  can  accept  a 
range  of  properties  as  input,  but  requires  only  clay  content  to  predict  porosity,  bulk  density,  elastic 
moduli,  and  compressional  and  shear  velocities.  The  model  has  been  validated  against  well  logs  from 
ODP  boreholes  on  the  Amazon  Fan. 

A  1-D  poro-elastic  model  of  seismic  wave  propagation  is  also  complete.  The  model  has  been 
exercised  in  simulating  synthetic  seismograms  of  experimental  strata  using  a  range  of  source 
frequencies  (e.g.,  100  -  2000  Hz)  and  source  signatures  (e.g.,  chirp,  airgun  and  watergun)  (Pratson  et 
al.,  2001)  (Fig.  1). 

Finally,  the  unconfined  debris  flow  experiments  needed  to  validate  the  2-D,  layer-averaged  debris  flow 
model  have  been  completed. 
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RESULTS 


•  The  shelf  model  predicts  that  continental  shelves  should  evolve  to  a  relatively  flat  gradient,  the  dip  of 
which,  if  any,  is  controlled  by  the  dissipation  of  wave  energy  across  the  shelf. 

•  The  slope  model  predicts  that  repeated  slope  failure  creates  an  abrupt  shelf  break,  while  slope 
bypassing  fosters  and  preserves  a  more  gradual  rollover  across  the  shelf  edge. 

•  Using  clay  content  as  the  only  input,  the  physical  property  model  can  predict  the  porosity,  bulk 
density  and  velocity  of  sediments  from  the  Amazon  Fan  to  within  a  RMS  error  of  <  10%. 

•  The  seismic  model  successfully  reproduces  the  greater  attenuation  of  seismic  energy  by  sandy 
sediments  versus  those  that  are  clay  rich  over  seismic  frequencies  commonly  used  to  image  shelf  and 
slope  strata  (Fig.  1). 

IMPACT/APPLICATIONS 

•  The  prediction  that  continental  shelves  should  evolve  to  relatively  flat  surfaces  is  due  to  the  shelf 
model's  simplicity,  but  the  result  raises  the  interesting  question  of  what  creates  and  controls  the  slope 
of  the  shelf  (e.g.,  modem  marine  processes  or  past  subaerial  processes  during  sea  level  lowstands 
like  rivers?).  This  is  a  question  that  Dr.  Pat  Wiberg  and  others  are  pursuing. 

•  Based  on  observations,  shelf-break  abruptness  appears  in  general  to  be  inversely  correlated  with 
sediment  supply.  Abrupt  shelf-breaks  are  characteristic  of  stationary  continental  margins  with 
moderate  to  low  sediment  supply,  while  rounded  shelf-breaks  are  common  along  prograding  margins 
with  high  sediment  supply  (O’Grady  et  al.,  2000). 

•  While  clay  content  can  be  used  to  predict  seismically  important  physical  properties  of  marine 
sediments,  it  itself  is  a  property  that  is  poorly  measured,  if  it  is  measured  at  all. 

•  With  some  pre-existing  knowledge  of  sediment  type,  it  appears  that  a  1-D  poro-elastic  model  can  be 
used  to  predict  what  seismic  source  would  yield  the  best  images  in  a  given  area.  Conversely,  under 
certain  conditions  (e.g.,  gas-free  sediments),  the  model  may  allow  for  gross  predictions  about 
sediment  type  from  the  attenuation  of  a  given  seismic  source. 

TRANSITIONS 

The  linked  coastal-plain,  shelf,  slope  and  rise  models  have  been  provided  to  Dr.  Mike  Steckler  (L- 
DEO).  He  has  already  incorporated  them  into  SEQUENCE,  and  he  has  used  them  to  produce  example 
simulations  of  margin  evolution  (Steckler  et  al.,  2001).  Similarly,  the  physical  property  and  seismic 
models  have  been  provided  to  Dr.  James  Syvitski.  He  will  be  using  them  in  the  future  to  generate 
synthetic  seismograms  of  strata  simulations  produced  by  SEDFLUX.  Dr.  Jim  Buttles  has  been 
working  with  the  results  from  the  debris  flow  experiments  in  testing  the  debris  flow  model.  And  the 
experiment  results  have  also  been  used  by  a  Duke  undergraduate  who  analyzed  them  for  his  honor 
thesis  this  past  spring  (2001)  (Henderson,  2001) 
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RELATED  PROJECTS 


I  am  continuing  my  work  with  Dr.  Dave  Cacchione  on  the  relationship  between  the  dynamics  of 

internal  waves  and  the  gradients  of  the  continental  slopes  (Cacchione  et  ah,  in  2001).  We  are  now 

extending  our  analyses  beyond  the  STRATAFORM  study  areas,  and  looking  at  the  relationship  on  a 

global  basis. 
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Figure  1.  Synthetic  seismograms  of  experimental  strata  containing  buried 
channels:  (top)  strata  are  modeled  as  being  predominantly  sandy;  (bottom)  strata 
are  modeled  as  being  predominantly  muddy.  Note  difference  in  attenuation. 
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